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ABSTRACT 

Geometrical optics provides an excellent description for quasar images crossing caus- 
tics which are formed by gravitational microlensing of objects like Q2237+0305. Within 
this approximation the source size can be estimated from the maximum magnification 
reached at caustic crossings. We evaluate the limitations imposed by diffraction on caus- 
tics using the formalism developed by Ulmer & Goodman (1995). Close to a caustic a new 
characteristic length, smaller that the Fresnel length, enters the problem, limiting the an- 
gular resolution to about 0.2 pico arcsecond, or equivalently ~ 3 x 10 9 cm at the source. To 
achieve this resolution the brightness must be monitored at time intervals of a few seconds. 
If a significant fraction of quasar luminosity comes from sources smaller than those limits 
then interference effects would make the observed intensity oscillate, in a close analogy 
with a two slit experiment. The characteristic period of such oscillations is expected to 
be about one tenth of a minute. If such oscillations are detected then photometry carried 
out at a single site may permit the determination of the caustic transverse velocity, and 
therefore may permit a direct conversion of the time units of brightness variations to the 
linear units at the source. 

Subject headings: Gravitational lensing - dark matter - quasars: structure -quasars: Q2237- 
+0305 



- 2 - 



1. INTRODUCTION 

The quadrupole gravitational lens 2237+0305 (Huchra et al. 1985) is the first, and 
so far the best studied case of gravitational microlensing at large optical depth. The four 
macroimages are formed by the whole galaxy acting as an astigmatic lens. The four lines 
of sight pass through four different inner regions of the lensing galaxy and are subject to 
microlensing by the stars located near those lines of sight. Some of the four macroimages 
were observed to vary on a time scales of months (Racine 1992 and references therein) 
and even days ( Pen et al. 1994). The optical depth to macrolensing combined with 
the shear make the magnification of all four macroimages so large (cf. Wambsganss & 
Paczyhski 1994, and references therein) that the microlensing is due to a collective effect of 
many stars, as emphasized by Wambsganss (1992). This is best seen at the magnification 
patterns, or equivalently the illumination patterns which clearly demonstrate that high 
magnification events are mostly caused by the source crossing one of the many caustics 
created by the microlensing, and almost none is caused by an isolated star (Wambsganss, 
Paczyhski & Schneider 1990). 

The generic structure of almost every caustic crossing is the following. On one side of 

1/2 

the caustic the total magnification of a point source varies as A = A min + (d c /r s ) 1 (cf. 
Fig. 1), where A min is the magnification on the other side of the caustic, r s is the linear 
distance between the source and the caustic, and d c is a characteristic length related to 
the caustic. In this paper all lengths, distances, and velocities in directions perpendicular 
to the line of sight are measured in the lens plane, while their equivalents in the source 
plane are denoted by the same symbols with hats. The length scale d c is proportional to 
the Einstein ring radius r E of individual micro lenses, the stars in the lensing galaxy. If the 
source moves with a velocity V at an angle 9 to the caustic then the magnification varies 
with time as 

A « A min + (dc/AtVfanO) 1 ' 2 = A min + (At c /At) 1/2 , (la) 

At c = d c /V sin 9 = ar E /V sin 9 ~ M 1 ^ 2 /V sin 9, a = d c /r E , (lb) 

where At c is the characteristic time related to the caustic, M is the mass of microlensing 
stars, and a is a dimensionless parameter of the order unity. The magnification of a point 
source increases all the way to infinity upon crossing any caustic. When the caustic is 
crossed the image magnification drops discontinuously to a finite value A m in- Naturally, 
the crossing is equally likely to be in the opposite direction, with a discontinuous increase 
of the image magnification from A min to infinity, followed by a gradual decline described 
with the equations (la,b). 

If the source has a finite size, d s , then the maximum magnification at the caustic 
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crossing is 



A 



-4min + (d c /d s ) 



1/2 



(2) 



max 



and the discontinuous jump in the magnification is replaced by a change on a time scale 



where At s is the characteristic time related to the source. 

Note that the high magnification events caused by caustic crossings are characterized 
by two distinctly different time scales. One time scale, At c is proportional to the square 
root of some average mass of the stars that are responsible for microlensing. The other 
time scale, At s is proportional to the source size. With frequent and accurate photometric 
measurements it should be possible to partly deconvolve the source structure (Grieger, 
Kayser, & Refsdal 1988; Wambsganss & Paczyhski 1991). 

Pen et al. (1994) presented the first observations of 2237+0305 variability which 
were interpreted as caused by a caustic crossing. The image A increased its brightness by 
~ 1.5 magnitude over two months, corresponding to At c »3x 10 6 s, and subsequently 
its brightness dropped by ~ 1.5 magnitude in less than a week, corresponding to At s xs 
3 x 10 5 s. If our transverse velocity with respect to the lensing galaxy at z g = 0.04 is 
~ 10 3 km s _1 , then the projected velocity at the source at z s = 1.7 is V ~ 10 4 km s _1 , 
and the source size may be estimated to be d s ~ V6At s ~3x 10 14 cm. This is just ~ 20 
astronomical units, much less than the estimates based on accretion disk models: d s ~ 
1.5 x 10 16 cm (Rauch & Blandford 1991) or d s ~ 4 x 10 15 cm (Jaroszyhski, Wambsganss, 
& Paczyhski 1992, cf. also Czerny, Jaroszyhski, & Czerny 1994; Jaroszyhski & Marck 
1994; Witt & Mao 1994). 

There is a possible way to reconcile the results of Pen et al. (1994) with the standard 
disk model of the quasar 2237+0305 if the event was caused by a very fast moving caustic 
crossing the image of the quasar. Such fast caustics were found by Kundic & Wambsganss 
(1993) and by Kundic, Witt & Chang (1993) to be the consequence of random motion of 
the stars responsible for microlensing. Only future observations and the well sampled light 
curves are likely to resolve this issue. 

In the geometrical optics approximation the caustics are infinitely sharp and in prin- 
ciple the resolution of quasar structure is limited only by the accuracy and frequency of 
photometric measurements during the rapid change of the apparent brightness caused by a 
caustic crossing. However, at some level diffraction effects must limit the resolving power 
of any microlensing system. The aim of this paper is to determine what is the highest 
resolution that might be achieved. In the next section we treat the problem theoretically 
and in Sec. 3 we apply our calculations to Q2237+0305. Discussion of the observability of 
the effects follows in the last section. 



At s fa d s /V sin#, 



(3) 
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2. DIFFRACTION BY CAUSTICS 

The diffraction effects near optical caustic have the same nature as in a two-slit ex- 
periment: they are caused by the interference of each photon passing the two paths corre- 
sponding to the two images which appear or disappear when the source crosses the caustic. 
The interference of radiation coming to an observer along different paths through the gravi- 
tational field, has been investigated by several authors. (Schneider & Schmidt-Burgk 1985; 
Deguchi & Watson 1986; Peterson & Falk 1991; Gould 1992; Stanek, Paczyhski & Good- 
man 1993; Ulmer & Goodman 1995). In most cases the calculations were limited to the 
simplest case of a single point-mass lens. The effects of physical optics in gravitational 
lensing are referred to as femtolensing (Gould 1992). 

Recently Ulmer & Goodman (1995, hereafter UG) have presented the method of 
calculating femtolensing effects for a general gravitational lens system. In particular they 
obtained the solution for the case of a source close to a caustic. In this paper we estimate 
the importance of the femtolensing for quasars. We are interested in the effects detectable 
with a broad band photometry, rather than the femtolensing effects in the spectra. 

We use a general expression for the time delay along a ray crossing the lens surface 
at point r = (x,y) and coming from a source point r s = (x a , y s ) as projected onto the lens 
plane (cf. Blandford & Narayan 1986): 

1 + z 

ct(x, y, x B , y a ) = ^ ((x - x s ) 2 + (y - y s ) 2 ) + (1 + z g )cr g (x, y) (4) 

where r is the total time delay and r g is the delay caused by the gravitational field of the 
lensing galaxy, D = D g D gs /D s (in standard notation) is the characteristic gravitational 
lens distance and z g is the redshift of the lens. Suppose we choose our coordinate system in 
such a way, that locally, near the origin, the line x = is tangent to a critical line related 
to a fold caustic. In the simplest case we have r xx = on the critical line x = 0. For 
a system consisting of point lenses one has r g)XX + r g ^ yy = everywhere, except singular 
points located at point masses. This equation fixes the value of the r yy derivative on the 
critical line, where r xx = 0. Introducing ip - a quantity proportional to the time delay and 
expanding it in the lowest nontrivial order we have: 

/ DcT loo 2 ^ 2 /k\ 

w = = -ax + y — x s x — y s y H — x„ H — y„ (5) 

v l + z g 6 2 s 2 ys K ' 

where a is defined by the third x-derivative of the time delay and can be made positive by 
the transformation x — > —x. There are other parametrizations possible near a critical line 
(i.e. Schneider & Weiss 1986; Kayser & Witt 1989; Witt 1990; Witt, Kayser, & Refsdal 
1993). The caustic line in our case is locally given by x s = 0. Using equation (la) we have 
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near the caustic: 

A-A^^^fJ , (6) 

where r s is a distance of the source from the caustic line and d c is a parameter characterizing 
the caustic. In our case r s = x s . The magnification due to the caustic related images is 
proportional to the inverse of the HVVjll determinant, where is given by equation (5). 
Comparing with the above formula we get: 

and substituting into the expression for ip we have: 

^ = 4&d£ X * + y2 ~ XsX ~ VsV + + 2 V *' ^ 

Physical optics calculations introduce the Fresnel length, which in our case can be 
defined as: 

f\D\ 1/2 fcD\ 1/2 

where A is the wavelength and u the frequency of light. With this definition one can 
calculate the complex amplitude (UG): 

* ( *"' Vs) = 2< dX dV 6XP U ) ' (10) 

The normalization is such, that in the absence of the lens ^(x s ,y a ) = 1. In the case of a 
point source near the caustic line we get: 

|vl>(x s ,y s )| = 2 1 /3 7r i/2 ^ /3 |Ai(-r B /df)| , (11) 

where Ai is the Airy function (Abramovitz & Stegun 1964) and the characteristic scale for 
its argument changes is: 

dp ( dp^ 1 " 



* = **{i) ■ (12) 

Both results are in agreement with UG but we use different parametrization. Using asymp- 
totic expansion of the Airy function one can see, that for r s » d{ the flux from a point 
source F(r s ) = |\1/| 2 — > (dc/rg) 1 ^ 2 as in geometrical optics calculation. 

We have checked the sensitivity of the interference effects to the size of the source and 
the width of the frequency filter used for observations. For a monochromatic incoherent 
source with surface brightness I(x s , y a ) one finds its flux by the convolution: 

F~ J dx s J dy s I(x s ,y s ) \^(x Sl y s )\ 2 (13) 
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In Figure 2 we show light curves of monochromatic sources crossing the caustic. All sources 
are of the gaussian shape with characteristic sizes of 0.1, 0.3, 1, 3 and 10 x df. For sources 
with sizes d s < d{ the interference patern is clearly visible and present at any distance 
from the caustic. The asymptotic expansion of Airy function (Abramovitz & Stegun 1964) 
Ai(— x) ~ sin(x 3 / 2 + 7r/4) where we neglected the slowly changing amplitude, shows, that 
the spatial distance between consecutive fringes is 

A, = f*g) I/a ("> 

1/2 

The spacing Ar s is proportional to the smoothed magnification, A ~ (d{/r a ) . For a 
source moving with the velocity V as measured in projection onto the lens plane and 
perpendicular to the caustic, the characteristic time for time variations is given as 

At EE ^ (15) 

This spatial and temporal characteristics of interference fringes is a generic feature of a 
motion across any caustic. 

Let us consider now a typical photometric measurement using a broad band filter with 
a band pass Auj/uj = 0.2 We assume the source spectrum to be flat in the filter range and 
use the "top hat" filter. We make convolution in frequency as well as in space. (Amplitude 
depends on frequency through dp and df). The results are shown on Figure 3 for the 
same shapes of sources as previously. Only close to the caustic (where the phase of \& is 
the same for all frequencies) the interference pattern can be noticed. Farther away the 
fringes disappear. The number of clear fringes is ~ uj/ Auj. 

3. APPLICATION 

The caustic patterns were modeled by Wambsganss (1990), Witt (1990), Wambsganss, 
Witt, & Schneider (1992), and Witt et al. (1993, hereafter WKR) among others. The last 
of the quoted papers can be applied directly to our problem, since it gives the useful 
characteristics of caustics in the lensing system of Q2237+0305. We use the image A 
characteristics in our parameter estimation. The dependence of source magnification near 
caustic on the average microlens mass given by WKR shows that the relevant scale is the 
Einstein radius calculated for the average microlens mass: 

r E = (^F^) V2 = 1.6 x 10 16 cm (M av /M ) 1 ^ 2 h^ 2 , (16) 

where /175 is the Hubble constant in units of 75 km s -1 Mpc -1 . Now using WKR Table 4 
and translating their notation to ours we get the value of a = (K 2 )/(m) where K is their 
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flux factor and (m) is the averaged microlens mass in solar units. For image A we get 
a = 0.8. The value of a may be somewhat different for other macro-image parameters. 
Substituting we have: 

d c = ar E = 1.28 x 10 16 cm (M av /M ) 1/2 /i~ 1/2 . (17) 

The Fresnel length for an optical wavelength calculated for the same system is 

d F = 6 x 10 10 cm A^o fc- 1/2 , (18) 

where A 50 oo is the wavelength in units of 5000 A. Using this estimates of the characteristic 
scales we find that the limiting size for a source to show the diffraction effects at optical 
wavelength is: 

I <di = ^dt^3x 10 9 cm (M av /M )" 1/6 \% 3 00 h^ 1/2 , (19) 

This is the upper limit to the source size as measured in the source plane, which for 
Q2237+0305 is at 8 times the distance to the lens plane. Also, this is the best resolution 
of a quasar that can be achieved through gravitational microlensing. 

A characteristic time of femtolensing variability (assuming the relative velocity of 
600km s -1 at the lens plane which translates to 5000km s _1 at the source plane, values 
used in Q2237+0305 modeling) is 

At » 6 s (M av /M )" 1/6 \% 3 00 h~ 5 1/2 . (20) 



4. DISCUSSION 



In our calculations we neglected the flux coming to an observer from images not related 
to the critical line. This is justified as the source crossing of a caustic does not influence 
other images and the pair of merging images dominates the total flux. 

The diffraction limits the resolution achievable through photometric monitoring of 

/S q 1/6 2/ 3 1/2 

caustic crossing to di w 3 x 10 cm (M av /M Q )~ 1 A 5 q 00 h 75 where M av is the average 
mass of microlensing objects in the galaxy that macrolenses Q2237+0305. This corresponds 
to the angular resolution of ~ df/D S: where D s is the angular diameter distance to the 
source, i.e. the resolving power of a caustic crossing in the Huchra's lens is ~ 10~ 18 radians 
or 0.2 pico arcseconds in optical light for stellar mass microlenses. This is a very impressive 
resolving power indeed. 

If there is a substructure in the Q2237+0305 with a scale smaller than df then inter- 
ference pattern should be detectable photometrically during the caustic crossing events: 
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a series of luminosity fluctuations with increasing or decreasing amplitude with the char- 
acteristic time scale of few seconds (cf. eq. [20]). Notice that if such fluctuations were 
detected they would allow the determination of the transverse velocity of the caustic, or 
equivalently they would allow the determination of a relation between the time scale of 
light variability and the linear size at the source using single site observation. If no inter- 
ference effects are detected the determination of the caustic velocity requires at least two 
observing sites separated by about one astronomical unit, as first pointed out by Grieger 
et al (1988). 

Our estimates demonstrate that in the case of Q2237+0305, and similarly other 
quasars, the geometrical optics is a very accurate approach. Investigating light curves 
of quasars crossing caustics one can probe their spatial structure down to scale as small as 
O.O5i?0. Currently we have no reason to expect that any quasar has structure as small as 
a tenth of the solar radius, but it is good to know that structure down to such a small scale 
can be resolved with gravitational microlensing. It is also important to realize that the 
variation of image brightness should be monitored on the shortest time scales allowed by 
current technology in order to determine observationally what is the highest magnification 
present at caustic crossings, and what is the shortest time scale on which the brightness 
varies. 

Adopting the geometric optics approximation we can estimate the source size d a com- 
bining the equations (lb) and (2): 

a 

d s = — — ~ 2 r E , (21) 

where a ps 1 (cf. Table 4 of WKR), te is the Einstein ring radius of the microlensing 
masses, A max is the maximum magnification reached at the caustic crossing, and A m ; n is 
the magnification just outside the caustic. Notice, that the estimate does not require any 
knowledge of either the transverse velocity V, or the angle 9 (cf. eq. [la]). 

We thank Hans Witt for the discussion. This project was supported with the NSF 
grant AST 93-13620, NASA grant NAG5-2759 and KBN grant 2-P03D-020-08. 
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FIGURE CAPTIONS 

Fig.l. A schematic representation of a generic light curve when a point source crosses a 
caustic, as expected within a geometric optics approximation. 

Fig. 2. Caustic-related magnification for gaussian shaped, monochromatic sources as a 
function of their distance from the caustic. The plots are for the sources with 
the characteristic size of 0.1 (the thinnest solid line), 0.3, 1, 3 and 10 x <i f (the 
thickest line), where df is characteristic length near the caustic (eq. [12]). The 
geometrical optics result for a point source (eq. [6]) is shown for comparison as 
a dotted line. We adopt parameters of Q2237+0305 lensing system (eqs. [17,18]) 
for M av = 1 M & and A 50 oo = 1 to get the magnification values near the caustic. 

Fig. 3. The same as on Fig. 2 but for sources seen through a broad filter Alu/uj = 0.2. 



